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1. FUNDAMENTAL AND DERIVED UNITS 


Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbreviation 

Length 

Time - _ 

Force. 

l 

t 

F 

meter 

second 

weight of 1 kilogram. _ 

m 

s 

kg 

foot (or mile) 

second (or hour) 

weight of 1 pound 

ft (or mi) 
sec (or hr) 
lb 

Power - — 
Speed.. _ __ 

P 

V 

horsepower (metric) 

f kilometers per hour 

\meters per second. _ 

kph 

mps 

horsepower 

miles per hour — — 

feet per second 

hp 

mph 

fps 


Weight =mg 

Standard acceleration of gravity =9. 80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass= — 


2. GENERAL SYMBOLS 

v 


Moment of inertia =mP. (Indicate axis 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


of 


Kinematic viscosity 
p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m -4 -s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft -4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


Area 

Area of wing 
Gap 
Span 
Chord 

. . 6 2 
Aspect ratio, -g 

True air speed 

Dynamic pressure, pV 2 

Lift, absolute coefficient Cl = jjg 

D 

Drag, absolute coefficient C D — ^g 
Profile drag, absolute coefficient Co^^g 
Induced drag, absolute coefficient C Di 


Parasite drag, absolute coefficient C 'o v — ^g 


3. AERODYNAMIC SYMBOLS 

i w Angle of setting of wings (relative to thrust line) 
i, Angle of stabilizer setting (relative to thrust 
line) 

Q Resultant moment 

Q Resultant angular velocity 

R Reynolds number, p — - where l is a linear dimen- 

sion (e.g., for an airfoil of 1.0 ft chord, 100 
mph, standard pressure at 15° C, the corre- 
sponding Reynolds number is 935,400; or for 
an airfoil of 1.0 m chord, 100 mps, the corre- 
sponding Reynolds number is 6,865,000) 
a Angle of attack 

c Angle of downwash 

a* Angle of attack, infinite aspect ratio 

a t Angle of attack, induced 

Dt oi a Angle of attack, absolute (measured from zero- 

<L$> lift position) 

y Flight-path angle 


Do 


C 


Cross-wind force, absolute coefficient Cc — ^ 


r 


V 
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TABLE II.- TEST DATA 


I Nozzle number 

Test gas 

H 

Jh 

d> 

l 

0 

1 

Test— chamber 
static pressure 
(microns Hg) 

Stream static 
pressure 
(microns Hg) 

Stagnation 
pressure 
(mm of Hg) 

Stagnation 

temperature 

(°p) 

Test— chamber 
surface temperature 

Model temperatures 


Iron end (°F) 

Center point 
(°F) 

Constantan end 
(°F) 

Model drag 
(milligrams ) 

1 

n 2 

2.49 

71.0 

72.8 

2.70 

73.0 

72.6 

101.3 

127.7 

95.0 

12.2-12.1 


n 2 

2.78 

120.0 

120.2 

5.00 

73.0 

72.2 

94.0 

119.0 

87.2 

21.7 


He 

• 55 

18.0 

19.5 

.70 

70.2 

70.0 

80.0 

87.0 

78.5 

1.0-0. 8 


He 

2.08 

51.0 

49.4 

2.15 

72.2 

70.2 

123.7 

I89.O 

117.0 

4. 8-4. 5 


He 

2-37 

89.0 

95.6 

4.00 

64.0 

62.2 

120.0 

211.0 

100.2 

15.5-16.0 

2 

n 2 

2.07 

71.6 

74.4 

1.20 

75.5 

75.2 

103.3 

131.7 

96.5 

8.8 


n 2 

2.25 

119.0 

132.7 

2.20 

74.7 

75.2 

95.7 

122.7 

90.2 

18.5-17.5 


He 

1.57 

53.5 

52.0 

1.00 

75.5 

75.5 

105.7 

141.5 

98.0 

4. 8-4. 7 


He 

1.75 

91.0 

91.0 

1.70 

65.0 

67.5 

114.4 

183.3 

100.6 

8.7 

3 

n 2 

1.85 

68.7 

68.7 

.90 

75.0 

74.4 

101.7 

129.5 

95.5 

6. 8-6. 7 


He 

.73 

32.0 

33.0 

.35 

76.2 

74.0 

82.5 

93.7 

78.5 

1.7-1. 5 


He 

1.00 

44.3 

46.5 

.60 

67.4 

68.7 

87.6 

110.5 

81.4 

2.7 

4 

He 

.89 

22.6 

22.6 

.19 

68.5 

74.7 

80.7 

87.5 

78.2 

1.1-1. 2 


He 

1.67 

55.0 

55.0 

.50 

74.0 

74.4 

101.0 

132.0 

93.7 

3. 5-3. 6 


He 

1.68 

78.5 

78.5 

.70 

68.2 

69.4 

107.0 

162.5 

95.5 

7.0-6. 7 


He 

1.73 

102.0 

102.1 

.90 

70.7 

70.8 

110.7 

174.5 

97.5 

10.4-10.0 


He 

1.79 

117.0 

122.1 

1.05 

71.2 

71.5 

108.7 

175.7 

95.0 

12.6-12.4 


He 

1.83 

130.0 

135.0 

1.15 

70.8 

71.2 

106.3 

173.5 

92.4 

14.7-14.5 

5 

He 

.78 

18.6 

18.6 

.15 

73.7 

74.3 

76.6 

80.5 

74.8 

0.75 


He 

1.23 

43.5 

46.0 

.30 

76.3 

74.3 

94.2 

117.7 

88.2 

2 . 85-2 . 7 


He 

1.36 

63.6 

62.5 

.40 

74.2 

74.3 

102.7 

140.3 

94.5 

4.4— 4.2 

. 1 

He 

1.44 

80.5 

80.5 

.47 

73.3 

74.8 

105.0 

147.8 

97.0 

6. 0-5. 8 


1 Average Mach number over center half inch of stream. 
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A COMPARISON OF THEORY AND EXPERIMENT FOR HIGH-SPEED FREE-MOLECULE FLOW 1 

By Jackson R. St alder, Glen Goodwin, and Marcus 0. Creager 


SUMMARY 

A comparison is made of free-molecule-flow theory with the 
results of wind-tunnel tests performed to determine the dray and 
temperature-rise characteristics of a transverse circular cylinder. 

The model consisted of a 0.0031 -inch-diameter cylinder 
constructed of butt-welded iron and constantan wires. Tempera- 
tures were measured at the center and both ends of the model. 
The cylinder was tested under conditions where the Knudsen 
number of the flow (ratio of mean-free-molecular path to 
cylinder radius ) varied between 4 and 185. The molecular 
speed ratio (ratio of stream speed to the most probable molecular 
speed ) varied from 0.5 to 2.3. In terms of conventional flow 
parameters , these test conditions correspond to a Mach number 
range from 0.55 to 2.75 and a Reynolds number range from 
0.005 to 0.90. 

The measured values of the cylinder center-point temperature 
confirmed the salient point of the heat-transfer analysis which 
was the prediction that an insulated cylinder would attain a 
temperature higher than the stagnation temperature of the 
stream . 

Good agreement was obtained between the theoretical and the 
experimental values for the drag coefficient. As predicted , 
theoretically , the drag coefficient was independent of Knudsen 
(or Reynolds number) variation and depended primarily upon 
the molecular-speed ratio. 

INTRODUCTION 

The various phenomena connected with the high-speed 
flows of rarefied gases have recently been the subject of 
several analytic investigations (references 1, 2, 3, and 4). 
Interest in this phase of aerodynamics and thermodynamics 
has been aroused primarily because of the necessity of 
being able to compute forces on and heat transfer to missiles 
which may fly at great altitudes at high speed. In another 
analytic exploration of this field, simplified concepts of 
aerodynamic and thermodynamic phenomena occurring 
at high altitudes have been used by Whipple (reference 5) 
in an investigation of the properties of the upper atmosphere 
hy means of measurements of the brightness and deceleration 
of meteors. 

Experimental work in the field of high-speed rarefied gas 
flows lias been practically nonexistent up to this point. The 
pioneer work of Epstein (reference 6) on the drag of spheres 


in a rarefied gas was concerned only with very low-speed 
motions. Consequently, the purpose of the present in- 
vestigation was twofold: (1) to obtain experimental data 
that could be compared with the predictions of an analytic 
treatment of aerodynamic and thermodynamic processes 
occurring at very high altitudes; and (2) as a means to this 
end, to develop a wind-tunnel and related instrumentation 
that would permit testing of models under simulated condi- 
tions of high-altitude high-speed flight. In connection 
with the second objective — the development of a wind 
tunnel and instrumentation — it should be pointed out that 
a low-density wind tunnel is being developed concurrently 
at the University of California and the authors wish to 
acknowledge the assistance they have received from the 
reports issued by this group, as well as the stimulus received 
from conversations with the University of California person- 
nel. 

At the beginning of the experimental work, it was decided 
that the initial model tests in the low-density wind tunnel 
should be confined to the free-molecule-flow regime despite 
the fact that this regime is probably of less technical im- 
portance than the slip-flow and intermediate regimes. This 
decision was reached because the free-molecule regime is 
susceptible to analysis; whereas the slip-flow and interme- 
diate regimes have not yet been adequately treated (cf., 
reference 2). In view of the many uncertainties connected 
with the development of entirely new types of wind-tunnel 
instrumentation, it was thought to be essential to test in a 
field where the experimental work could be guided by a 
theoretical background. 

The selection of the test body was governed both by the 
requirement that free-molecule flow be obtained and by the 
practical considerations arising from tunnel size and instru- 
mentation. The former condition required that the signifi- 
cant body dimension be small compared to the mean-free- 
molecular path. 

A cylinder with axis perpendicular to the stream was 
chosen as the body which most nearly fulfilled the require- 
ments of the test. The radius of the cylinder could be made 
small compared with the mean-free path. In addition, 
since it was desired to measure both aerodynamic forces and 
body temperature, the cylinder had the virtue of being 
doubly symmetrical so that only drag force had to be 
measured. The problem of temperature measurement could 


• Supersedes NACA TN 2244, “A Comparison of Theory and Experiment for High-Speed Free-Molecule Flow” by Jackson R. Stalder, Glen Goodwin, and Marcus O. Creager, 1950. 
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be solved by constructing the model from dissimilar metals 
so that the model itself formed a thermocouple. 

The present investigation is a first step toward an attempt 
to gain an understanding of the comparatively new and 
unexplored field of high-speed low-density flows. The ex- 
periments reported herein were made to provide a comparison 
with some important results of free-molccule-flow theory 
and to facilitate the solution of problems connected with the 
design, operation, and instrumentation of a low-density 
wind tunnel. The tests were conducted in the Ames low- 
density wind tunnel. 

NOTATION 

A surface area, square feet 

Cr I 

cj. components of total velocity, feet per second 

?z j 

C D total drag coefficient, dimensionless 

C D impinging drag coefficient , dimensionless 

C D re-emission drag coefficient, dimensionless 

D total drag force, pounds 

Ei rate of incident molecular energy, foot-pounds per 

second 

E r rate of re-emitted molecular energy, foot-pounds per 
second 

E r rotational energy, foot-pounds per molecule 
E t rate of incident translational molecular energy on 

front side of body, foot-pounds per second 
/($) dimensionless function of s defined by /(s) — Z] 

(s 2 +3)+Z 2 (V+2) 

f°(s) dimensionless function of s defined by f°(s)=Zi 

(s*- + 2)+Z 2 (s 2 +l) 

g(s) dimensionless function of s defined by </($) = 3 

(Zi +Z 2 ) 

g°(s) dimensionless function of s defined by g°(s) = 2 
(Z 1 +Z 2 ) 

G normal force, pounds 

7 0 modified Bessel function of first kind and zero order 

/j modified Bessel function of first kind and first order 

J mechanical equivalent of heat, 778 foot-pounds per 

Btu 

k Boltzmann constant, 5.66X10 -24 foot-pounds per 

°F per molecule 

K Knudsen number (rat io of mean-free-molecular path 

to cylinder radius), dimensionless 

H 

l v \ direction cosines, dimensionless 

u I 

L length, feet 

m mass of one molecule, slugs 

M Mach number, dimensionless 

n number of molecules striking body, molecules per 

square foot, second 

N number of molecules per unit volume of gas 


j> free-stream static pressure, pounds per square foot 

r cylinder radius, feet 

E universal gas constant , 1544 foot-pounds per pound- 

mole, °F absolute 

,<? molecular-speed ratio (ratio of stream mass velocity 

to most probable molecular speed), dimensionless 
T free-stream static temperature, °F absolute 

T 0 tunnel stagnation temperature, °F absolute 

T c cylinder temperature, °F absolute 

T u surface temperature, °F absolute 

T t test-section wall temperature, °F absolute 

l ' mass velocity, feet per second 

rJ 

f ' ;/ i components of mass velocity, feet per second 

/'J 

v r re-emission velocity, feet per second 

r m most probable molecular speed, feet per second 

11 molecular weight, pounds per pound-mole 

x\ 

y j Cartesian coordinates 

Zi dimensionless function of s defined by Zi = 7rc~ s2/2 

Io (S 2 / 2) 

Z 2 dimensionless function of s defined by 

Z 2 =7rsV’ 2 /2 l/ 0 ($ 2 /2) + 7 i (s 2 /2)] 

a accommodation coefficient, dimensionless 

ft reciprocal of most probable molecular speed, sec- 

onds per foot 
T gamma function 

7 ratio of specific heats, dimensionless 

e emissivity, dimensionless 

6 angle of body element with respect to free stream, 

radians 

p density of free stream, slugs per cubic foot 

a Stefan-Boltzmann constant, 3.74X 10 -10 foot-pounds 

per square foot, °F 4 

\jz dimensionless group defined by 

sin d [X+erj (s sin 0)] e s2 sin2 d 

1 + 

1 H“ V tt s sin 0 [1 + er/(s sin 0)] e s2 *™ 2 o 
\p f dimensionless group defined by 

1 — ' ^ ^ s sin 0 [1 — erf (s sin 0)] e 82 sin2ff 
1 — A tt .s‘ sin 0 [I —erf (s sin 0)] 

X dimensionless group defined by 

X = e-« 2 sin 2 0 ~j- ^/tt s sin 0 [1-f erf (.v sin 0)] 

X dimensionless group defined by 

X = f~* 2 sin 2 0 — ^ tt ,s* sin 0 [1 — erf (.*? sin 0)] 

erj(a ) error function e~ z2 dx^ dimensionless 
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ANALYSIS 

The general problem of the calculation of aerodynamic 
forces and heat transfer to a body in a free-molecule-flow 
field is amenable to solution by application of the simpler 
concepts of the kinetic theory of gases. The simplification 
arises from the assumption that the motions of molecules 
that strike the body are in no way affected by collisions with 
molecules that have struck the body and have returned to 
the main stream of gas. This basic postulate allows the 
molecular motions of the oncoming stream to be treated as 
having the classical Maxwellian velocity distribution super- 
imposed on a uniform mass velocity. Hence, the severe 
complication of dealing with a nonuniform-gas is avoided. 
Of course, this assumption does not imply that returning 
molecules do not collide with oncoming molecules since 
collisions do occur, on the average, at a distance from the 
body equal to one mean free path. However, if this mean 
free-path length is large compared with the body dimensions, 
it can be seen from geometrical considerations that the 
probability is small that the impinging molecules that do 
collide with returning molecules will strike the body without 
first having had their original Maxwellian velocity restored 
by collisions with other impinging molecules. It would 
seem probable, also, that a high value of mass velocity com- 
pared with the random thermal velocity of the gas would 
favor free-molecule flow because the re-emitted molecules 
would tend to be swept along with the stream away from the 
zone of influence of the body. 

In general, the mechanism of energy and momentum 
transfer for free-molecule flow are inseparable. In effect, 
this means that the drag on a body is affected by the surface 
temperature of the body. This state of affairs arises because 
of the fact that a portion of the momentum imparted to the 
body is due to the emission of molecules from the surface 
I and the velocity of emission depends upon the temperature 
of the body. This phenomenon will be discussed in detail 
later. 

Because of the assumption that the motions of the oncom- 
ing and re-emitted molecules can be treated independently, 
the computation of the drag on and heat transfer to a 
cylinder can be broken down into separate parts — the part 
due to impact by impinging molecules and a part due to 
! re-emission of the molecules from the surface. By the use 
of this basic principle, it is possible to calculate the drag 
and heat-transfer characteristics of a transverse cylinder 
in a free-molecule-flow field. 

HEAT TRANSFER TO A TRANSVERSE CYLINDER IN A FREE-MOLECULE- 

FLOW FIELD 

The calculation of the heat-transfer or energy-exchange 
process between the cylinder and the gas is made by using 
i the basic method described in reference 4. The general equa- 
tions of reference 4 are expressed in Cartesian coordinate 
form and can be summarized as follows: 

The translational molecular energy incident upon the 
front side of an elemental plane area dA inclined at an angle 
6 with respect to the stream mass velocity V is, for a 
monatomic gas, 


dE t = n 



( 1 ) 


where n, the number of molecules striking unit surface area 
per second, is given by 



(2) 


For the rear side of the elemental area, similar expressions 
apply. 

, / mU 2 


dE 




2 
Nv 




m f 

X 


(3) 

(4) 


The total molecular energy incident upon the front and 
rear sides of dA is, for a monatomic gas, given by equations 
(1) and (3), respectively, since only translational energies 
are involved. However, for a diatomic gas, the expression 
for the incident energy contains an additional term due to 
the rotational component of the molecular internal energy, 
and this component must be added to the translational com- 
ponent. The rotational energy per molecule of a gas com- 
posed of rigid dumbbell molecules is 


E R =kT (5) 

Then, the total incident molecular energy on tiie front 
side of dA is, for a diatomic gas, given by the following 
expression: 

dE t =dE t +n E R dA=n p ~ '+(*+ 1)4 7’J dA (6) 

and, for the rear side of dA , 

dE/=dE/+n'E R dA=n' pp-(*'+ l)*5rj dA (7) 

The method of calculation of the energy transported from 
the elemental area by re-emitted molecules is shown in ap- 
pendix A. The result derived therein can be applied to the 
calculation of the energy transport by a monatomic gas from 
the front side of the area by combining equations (Al) and 
(A 12) to yield 

dE r =dE t ( 1 — a) + 2 ctnkTcdA (8) 

and, for the energy transport by a diatomic gas, as 

dE r =(dE t -\-nk Fd H)(l — a)-\-3atikT c dA (9) 

Parallel equations can be written for the rear side of the 
area, dA. 

The net radiant energy exchange is (assuming gray-body 
radiation) 

e<r(T c *-T*)dA (10) 


where T t is the effective temperature of the surroundings to 
which the area dA is radiating. This expression implies 
that the body of which dA is a part is small in comparison 
with the surrounding volume. 

The problem now becomes one of applying the foregoing 
equations, which were derived for a plane area placed at an 
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arbitrary angle with respect to the stream, to a specific 
body — in this case a cylinder with axis perpendicular to the 
stream velocity. For a cylinder of length L , the elemental 
area mav be written as 

dA—rLdO (11) 


Then, if internal energy input to the cylinder per unit area 
is denoted as Q (in the usual experimental arrangement this 
would be electrical energy), an equation can be written 
which expresses an energy balance on the cylinder. For a 
monatomic gas, the energy-balance equation is 


2arL ^ j n ( +#r) dd+ J * ’ n' (~- + ^kT^ dd^-4arLkT c ( £ nd6+ £ n'<!6)+2irrL |^Q-e<7 ( 27— iy)J 

For a diatomic gas, a similar equation applies 


= 0 
( 12 ) 


| //p— ^ +($+l)kT^dd+ | ‘ + — SatrLkTc^ ndd+ £ n'dd^+2TnL^Q— — 


0 


The reduction of the integrals appearing in equations (12) 
and (13) is a straightforward but lengthy procedure. The 
methods used to evaluate these integrals are given in ap- 
pendix B. The final results are written here, for a monatomic 
gas, as 

9 7 J ( 7 . + Z 2 )-[z, (.s- 2 + 2)+ Z 3 ( «*+§)]+ 


T 


b<r(7’/-7V)-y]=l 

pv m a L J 


(14) 


and, for a diatomic gas, as 

:s If (Z, + Z 2 )-[z,(.7+ 3)+ )]+ 

o_3/2 r~ “I 

— e«(TS-T t *)-Q =( 

VVm<X L J 


= o 


(15) 


In order to facilitate computation, values of the dimensionless 
groups, 


and 


f ( s ) = Zi(s 2 + 3 ) + Z 2 (s 2 + ^ ^ > 
5(s)=3(Z 1 +Z 2 ), 

/»(«)= Z,(«*+2) + £(«*+§)» 


^°G 9 )— 2(Zi + Z 2 ), 
appearing in equations (14) and (15), are tabulated in the 
table of appendix B. 

Equations (14) and (15) then completely describe energy 
transfer occurring between a transverse cylinder and its 
surroundings in terms of the cylinder temperature, static gas 
t cm perat u re and pressure , molecular-speed ratio , accom- 
modation coefficient, the effective temperature of the sur- 
roundings to which the cylinder is emitting radiation, the 
emissivity of the cylinder surface, and the internal energy 
input to the cylinder. The equations are subject to the 
following restrictions: 

1. The surface temperature of the cylinder is assumed to 
be constant, circumferentially. 

2. The cylinder receives and emits radiant energy as a 
gray body; that is, the emissivity is independent of wave 
length. Also, the cylinder is assumed to be small compared 
to the surrounding volume. 

3. The internal energ} r of the gas stream in which the 
cylinder is located has the value (3/2) kT for a monatomic 
gas and (5/2) Z'Tfor a diatomic gas. 


(13) 


DRAG ON A TRANSVERSE CYUNDLER IN A FBEE-MOLECTJLE-FLOW FIELD 


The general problem of aerodynamic forces acting on 
bodies in free-moleeule flow has been treated by several 
authors (references 1, 2, and 3). The basic method of 
calculation was identical in each case. The method con- 
sists of the assumption of a gas with a Maxwellian velocity 
distribution superimposed upon an arbitrary mass velocity. 
The number of molecules having specified velocities that 
impinge upon an elemental area inclined at an arbitrary 
angle to the stream mass velocity is computed and, from 
this calculation, the force due to the transmitted momentum 
is obtained. To this momentum due to molecular impinge- 
ment, however, must be added the momentum imparted to 
the surface by molecules emitted from the surface. The 
magnitude and direction of the force imparted by emitted 
molecules depends, among other things, on the type of 
emission that is assumed to occur. A condensed discussion 
of the several possible types of surface emission can be found 
in reference 4. An extensive treatment of the subject is 
given in reference 7. 

In order to determine the total aerodynamic force on a 
body of arbitrary shape, the differential momentum imparted 
to an element of the body surface is computed, as outlined 
in the preceding paragraph, and the total momentum 
obtained by integration over the body surface. Also, since 
the impinging and re-emitted momenta can be added 
linearly, the resulting total momentum can be calculated 
for any desired type of re-emission. 

In this paper, completely diffuse scattering of the incident 
molecular stream from an insulated surface is assumed. 
By this assumption it is meant that (a) the molecules that 
are scattered from the surface are emitted with a Maxwellian 
distribution of speed, (b) the direction of molecular emission 
from the surface is controlled by the Knudsen cosine law 
(cf., references 7 and 8), (c) the body temperature is assumed 
to be uniform over the surface, and further (d) that con- 
servation of number of molecules is maintained for each 
element of the body surface. These assumptions seem to 
be in accord with physical fact, as discussed in reference 6. 

The specific problem of the aerodynamic drag on a trans- 
verse cylinder has been treated in references 1 and 3. The 
reason for repetition of the calculation here is that the 
results have been obtained in a compact closed form which 
is amenable to easy calculation and, further, that a solution 
for the re-emission drag is obtained which involves measur- 
able physical quantities. Appendices A and C contain the 
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details of the calculations and only the final results are 
presented here. The results of the computations are as 
follows: The total drag acting on an insulated transverse 
circular cylinder of radius r and length L is given by 


D=^(2rL)(C Di +C Dr ) (16) 

where C Di and C Dr are the drag coefficients associated with 
impinging and re-emitted momenta. The value of C D . is 
shown in appendix C to have the value 


C D- ' 


ir 1/2 e- * 2 '2 
s 


{ I 0 (s 72)+(i±j^) [/ 0 (s72)+/,(s72)J j 


(17) 


and the value of C Dr for an insulated cylinder is shown in 
appendix A to be, for a monatomic gas, 


Zi(s 2 +2)+Z 2 
2(Zx-j- Z 2 ) 

and, for a diatomic gas, 

r _ ^ -X l Z ' (S 2 + 3}+Zi ( <s2+ 1 ) 

D '~ 4.7 V ' 3 (Z t +Z 2 ) 

APPARATUS AND TEST METHODS 



(18) 


(19) 


THE W INI) TUNNEL 

The Ames low-density wind tunnel, in which the present 
tests were conducted, is an open-jet nonreturn-type tunnel. 
A sketch showing the major components of the wind tunnel 
is shown in figure 1. The test gas, from a bottled source, is 
throttled into a settling chamber and then passes through an 
axially symmetric nozzle which discharges a free jet into an 
18-inch cubical test chamber. The test chamber is connected 
through a surge chamber to four booster-type oil-diffusion 
I pumps which continuously evacuate the system. The 
diffusion pumps are connected in parallel and are backed by a 
! mechanical vacuum pump which discharges to the atmos- 
phere. Valves are provided at the inlet to each diffusion 
pump in order to allow the pumping system to be isolated 
when the test chamber is opened to the atmosphere. 

A traversing mechanism is installed in the test chamber 
which permits three degrees of motion relative to the gas 
stream. The mechanism is used to position models and 
instrumentation in the desired location within the gas stream. 

The nozzles used in the tests were axially symmetric with 
an outlet diameter of 2 inches and were machined from clear 
plastic. The nozzle contours were determined by the method 
outlined in reference 9. No boundary-layer correction was 
applied to the nozzle contours. 

After the system was evacuated, the flow rate of the test 
gas into the settling chamber was adjusted until the stagna- 
tion pressure, as measured by a McLeod gage, reached the 
desired value. 

The static pressure in the test chamber was read on another 
McLeod gage and, if the resulting pressure ratio was that 


required to establish flow at the particular Mach number 
desired for the run, the tunnel was considered ready for 
tests. In order to match tunnel conditions of static pressure 
and pressure ratio with those that existed at an earlier time, 
it was sometimes necessary to alter the amount of throttling 
at the pump inlets due to a slight variation of pumping speed 
from day to day. 

The pressure level in the tunnel was varied by increasing 
or decreasing the flow of the test gas into the system. The 
static and stagnation pressures were again checked to de- 
termine if sufficient pressure ratio was available to establish 
flow at the desired Mach number. This procedure was some- 
what complicated by the fact that at the lower pressure levels 
the boundary layer in the nozzle would tend to thicken, thus 
causing a decrease in the Mach number compared with that 
obtained at the higher pressure levels. 

The wind tunnel could be operated through a range of test- 
section pressures of from 20 to 135 microns of mercury for a 
given nozzle and over a range of Mach numbers from 0.5 to 
2.75 by using several nozzles. 

INSTRUMENTATION ANI) THE TEST MODEL 

The static pressures which existed in the wind tunnel were, 
in general, too low to be measured with a conventional 
manometer or Bourdon-type pressure gage; therefore, two 
types of gages widely used in high-vacuum work were 
utilized, namely, the McLeod and Pirani gages. For a 
complete description of these gages the reader is referred to 
any standard text on high-vacuum technique (e. g., reference 
10); however, some of the special features of the gages used 
will be described. 

The McLeod gage used as a primary standard in these 
tests had an initial volume of 146 milliliters connected to a 
closed tube of 1 .95-millimeter diameter. This arrangement 
gave a scale length of 51.4 centimeters for a pressure range 
of 0 to 3000 microns. (For the sake of brevity, the term 
micron is used throughout as a unit of pressure rather than 
the term microns of mercury.) The diameter of the closed 
tube was measured by the mercury pellet method to within 
±0.01 millimeter and the difference in height of the mercury 
columns could be read to the nearest 0.1 millimeter by use 
of a micrometer scale. At a pressure of 100 microns, t his 
gage could be read to within 1 percent with ease. 

The gage was connected to the test chamber with an 
18-millimeter-diameter glass tube to insure a short time lag 
between equalization of the system and the gage pressures. 
A dry-ice and acetone cold trap was provided between the 
gage and the test chamber to prevent condensable vapors 
from reaching the gage. 

In general, it is difficult to obtain a completely independent 
method of checking the accuracy of pressure measurements 
made with a McLeod gage for pressures below about 1 
millimeter of mercury. Above 1 millimeter of mercury a 
U-tube manometer, filled with a very low-vapor-pressure 
oil and having a high-vacuum system (of the order of 10 -5 mm 
of Hg) connected to one leg, may be used. The McLeod 
gage used in the tests was checked in this manner and the 
agreement between the oil manometer and the McLeod gage 
was within ±3 percent. 
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Figure 1.— Ames low-density wind tunnel. 


An indication of the accuracy of the McLeod gage used in 
the tests in the pressure range below 200 microns was obtained 
by comparing the pressures indicated by the test McLeod gage 
with those obtained with a laboratory-type McLeod gage 
having a range of 0 to 200 microns, which has been calibrated 
against a laboratory standard McLeod gage at the manufac- 
turer’s plant. The agreement between the test gage and 
the laboratory-type gage was within ±3 percent at a pressure 
of 100 microns. 

The Pirani gage used in the tests consisted of an open 
glass tube, which contained a heated filament. The gas 
pressure in the tube determined the temperature and thus 
the resistance of the filament. A temperature-compensating 
tube of identical construction was sealed off at a very low 
pressure. The measuring and compensating tube filaments 
formed two legs of a Wheatstone bridge; and the unbalance 
of the Wheatstone bridge circuit, as indicated by a self- 
balancing potentiometer, was a measure of the pressure in 
the open leg of the gage. The Pirani gage was calibrated 
against the 0- to 3000-micron-range McLeod gage as a 
standard . 

A drift of the calibration of the Pirani gage with time was 
noted and for this reason the Pirani gage was calibrated 


before and after every test run. A typical calibration 
curve is shown in figure 2. 

Examination of figure 2 reveals that the sensitivity of the 
Pirani gage used was nearly constant from a pressure of 10 
microns to a pressure of 80 microns. Above about 200 
microns this Pirani gage becomes too insensitive to be 
usable as a pressure-measuring instrument. 



Figure 2.— Typical calibration of Pirani gage in nitrogen gas. 
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Figure 3. -Impact-pressure-measuring system. 


Whenever a pressure-measuring device is connected to 
a vacuum system by a tube or small orifice, two troublesome 
effects are present — time lag and gassing. The time lag of 
the pressure gage is due to the fact that, as the pressure in the 
main part of the system is changed, a volume of gas must 
How out of or into the pressure gage and its connecting 
tubing. At the low pressures encountered during the tests, 
this time lag may be of the order of hours if the gage volume 
and connecting tubing are not properly proportioned. 
Gassing, that is, gas being given up or adsorbed by the sur- 
face of the pressure gage and its connecting tubing, has the 
effect of causing the pressure gage to indicate a different 
pressure than that which exists in the main part of the 
system. Outgassing effects can be reduced materially by 
proper cleaning of the surface and by heating the surface 
under a vacuum to drive off water vapor and other materials 
which cling to ordinary surfaces. 

Schaaf, Mann, and Cyr (references 11 and 12) have 
investigated this problem in considerable detail, and show 
that, if the assumption is made that gassing is a function of 
surface area only, there exists an optimum geometry of gage 
volume and connecting tubing which will result in minimum 
time lags. References 11 and 12 were used as guides in the 
design of the impact and static-pressure-measuring devices 
which were used in the tunnel calibration tests, and, in 
addition, the time lag of these pressure-measuring systems 
was measured. 

The impact pressures encountered during the wind-tunnel- 
calibration tests were too high to be measured with a Pirani 
gage and, because of the difficulty of obtaining a suitable 
flexible tube which would have a small time lag, the McLeod 
gage could not be used. These impact pressures, therefore, 
were measured with a glass U-tube manometer filled with a 
Iow-vapor-pressure oil. The manometer was located inside 

9G5027 — 52 2 


the test chamber with one leg of the manometer open and the 
other end fastened with a glass connection to a Js-inch- 
diameter spherical-head impact tube. Figure 3 shows a 
diagrammatic sketch of this arrangement. The difference 
in oil level between the two legs, as read by a coordinate-type 
cathetometer located outside the tunnel, was the difference 
between the test-chamber static pressure and the impact 
pressure. The test-chamber static pressure was measured 
by the 0- to 3000-micron McLeod gage described earlier. 

The time lag of the system shown in figure 3 was measured 
by setting the pressure in the test chamber at approximately 
3 millimeters of mercury and then evacuating the chamber to 
approximately 10 microns as quickly as possible. The time 
required for the oil level in the two legs to equalize was 
approximately 15 seconds. This indicates a very short time 
lag for the practical case as the rate of pressure change 
encountered in the wind-tunnel calibration was small com- 
pared to that measured during the time-lag test. The 
procedure used in the wind-tunnel-calibration tests was to 
continue to take readings of the oil level until steady state 
was established. 

The static pressure of the gas stream from the nozzle was 
measured with a 0.10-incli-diameter ogival probe con- 
nected directly to a Pirani gage as shown in figure 4. 
Three 0.0225-inch-diameter holes spaced 120° apart were 
located at a point on the ogive where the pressure is equal 
to free-stream pressure as calculated by the method given 
in reference 13. The position along the probe where the 
surface pressure is equal to free-stream pressure varies 
somewhat with Mach number; however, this variation is 
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not large and static pressures measured by this probe agreed 
with pressures measured with a 4° cone. 

No boundary-layer correction was applied to the ogive, 
or cone contour; consequently, it is possible that some error 
in the measurement of static pressure existed. Figure 5 
shows the main dimensions of the static-pressure probe used. 

In order that the Pirani gage shown in figure 4 could be 
calibrated against the McLeod gage, which was connected 
to the test chamber, a 4-millimeter-bore stopcock was glass- 
welded to the Pirani gage to allow calibration with a rel- 
atively large opening and attendant short time lag. As the 
stream static pressures were measured with the stopcock 
closed, the following tests were conducted to measure the 
order of magnitude of the time lag of the system. The 
test chamber was pressurized to 165 microns and suddenly 
evacuated to approximately 5 microns. The rate of pressure 
change as measured by the Pirani gage with the stopcock 
open and with the stopcock closed is shown in figure 6. The 
procedure was repeated with a rising pressure and the results 
are also shown in figure 6. 

It can be seen from figure 6 that the two curves, repre- 
senting the runs with the stopcock open and closed, are 
essentially congruent after 70 seconds for decreasing pressure 
and after 50 seconds for the increasing pressure runs. With 
the stopcock open, the time lag of the Pirani gage should be 
very small (references 11 and 12) and the curve of pressure 
as a function of time includes the time required for the 
pumping equipment to evacuate the test chamber, the time 
lag in the self-balancing potentiometer, and the thermal 
lag in the Pirani filaments. The curve of pressure as a 
function of time for the stopcock closed includes all these 
factors plus the additional time lag imposed by the small 
holes in the static probe. The pressure variation during the 
time-lag tests was much more severe than any encountered 
during the static survey tests; therefore it can be concluded 
that, by waiting for at least 2 minutes before recording data, 
all errors due to time lag would be eliminated. 

The test model. — The model tested was a 0.0031-inch- 
diameter right -circular cylinder fabricated from iron and 
constantan wires and had an exposed length of l % 2 inch. 

Note.— A ll dimensions given in inches. 


X 

( t 

0 

0 

0. 05 

0.019 

. 10 

.030 

. 15 

.051 

. 20 

.064 

.25 

.075 

.30 

.084 

.35 

. 091 

.40 

. 096 

.45 

. 099 

..50 

. 100 


« 0.4200 > 

I.c 

i 

5 500 

| O.IOOO d 

0.0465 d 

T ' 


3 

X 5> 

\\ r - 

'> u -3 holes 0.0225 diam 

\ - 1 


^ a 

Figure 5.— Static-pressure probe used in the tests. 



Figure G.— Time response of the static-pressure-measuring system. 

The iron and constantan wires were electrically butt-welded 
together to form a thermocouple junction. Figure 7 shows a 
microphotograph of the model. It can be seen that the 
welding operation did not appreciably alter the shape of 
the model. The butt-welded junction was located in the 
center of the model and a 0.0031 -inch iron wire was silver 
soldered to the constantan side of the model at a distance of 
r / l6 inch from the center junction. A 0.0031-inch-diameter 
constantan wire was similarly fastened to the iron side of 
the model. These two auxiliary thermocouple junctions 
were used to measure the temperature at the ends of the 
model. 

The entire model was covered with fine-grain soot ob- 
tained from an acetylene flame in order that the emissivitv 
would be constant over the total model length and so that 
an established value of emissivity of 0.95 could be used to 
evaluate the heat loss by radiation. 

Thermocouple voltages from the center and end thermo- 
couples were read on a null-type laboratory potentiometer 
having a least count of 0.01 millivolt. 

The drag balance. — A specially constructed microbalance 
shown in figures 8 (a) and 8 (b) was used to measure the 
drag of the model tested. The balance consisted of a beam 
mounted upon an agate knife-edge which rested in a steel 
V -block, as shown diagrammatically in figure 9. The model 
was mounted between two arms extending down from either 
side of the beam. Electrical leads were transferred from 
the moving beam to the balance frame through mercury 
pools. 

A drag force on the model caused an unbalance of the 
beam which was counteracted by a magnetic-restoring force 
produced by an armature located in a magnetic field pro- 
duced by an alternating-current coil. An optical lever 
approximately 2 feet long was used to determine the zero 
position of the balance. 
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Figure 7. — Micro-photograph of model showing weld at junction of iron and 
constantn n. Enlarged 00 times. 


The magnitude of the drag force in terms of magnetizing 
current was determined by hanging known weights on a 
lever arm attached to the balance beam and measuring the 
current necessary to restore the balance to zero position. 
Figure 10 shows a typical calibration of the drag balance. 

The model was mounted between the extension arms of 
the drag balance and was shielded over all but the center 
l %2 inch, as shown in figures 8 (a) and 8 (b). A check was 
made on the effectiveness of the shielding on the supporting 
members. A small shield was mounted to the center main 
shield (extension arms) covering the center portion of the 
model (normally exposed). The completely shielded model 
was then placed in the gas stream. No drag forces were 
detected. 

It was found that the balance exhibited an erratic calibra- 
tion shift from day to day, which necessitated calibration of 
the balance before and after each drag determination. 

CALIBRATION OF THE WIND TUNNEL 

The purpose of the calibration of the wind tunnel was to 
determine the Mach number and static-pressure distribu- 
tions of the stream in the plane at which the model was 



(a) Front view. 


Figure 8.— Photograph of drag balance showing force shields. 


(b) Rear view. 
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Figure 9.— Diagrammatic sketch of drag balance. 


located — / inch downstream of the exit plane of the various 
nozzles used. 

The method used was to survey the stream with the impact- 
and static-pressure probes described earlier. These impact- 
and static-pressure probes, together with their pressure- 
measuring gages, were mounted 2 inches apart upon a 
common support from the traversing mechanism and read- 
ings were taken at K 6-inch intervals from the center to the 
edge of the stream. This method of mounting allowed both 
the impact and static pressures to he measured during the 
course of one run, thus obviating the necessity of matching 
tunnel conditions on separate runs. 

Table' 1 lists the nozzles tested, their design Mach number, 
and the tunnel conditions for which the surveys were made. 

The static pressure varied across the stream issuing from 
each of the nozzles tested. For most of the nozzles, the 
pressure was somewhat higher near the center portion of 
the stream than was the pressure of the quiescent gas in 
the test chamber. Certain of these nozzles, however, did 


TABLE I. — WIND-TUNNEL-CALIBRATION DATA 


Nozzle 

number 

Design 

Mach 

number 

Test 

gas 

Test cham- 
ber static 
pressure (ir i- 
erons Hg) 

Stagna- 
tion pres- 
sure 

(mm Hg) 

Test ' 
Mach 
number 

1 

4. 5 

Na 

71.0 

2/7 

2. 49 



N 2 

121.0 

5. 0 

2. 78 



He 

18. 7 

.70 

.55 



He 

51.5 

2.2 

2. 08 



He 

90.0 

3.8 

2.37 

2 

3.6 

X 2 

70. 0 

1.2 

2. 07 



n 2 

117. 0 

2.2 

2. 25 



He 

51.7 

.95 

1.57 



He 

89.2 

1. 7 

1.75 

3 

3.0 

N 2 

71.0 

0. 90 

1.85 



He 

31.0 

.35 

.73 



He 

45. 0 

.55 

1.00 

4 

3. 11 

He 

20.0 

0. 19 

.89 



He 

54. 7 

.55 

1.67 



He 

79.5 

. 70 

1.68 



He 

103.0 

.90 

1.73 



He 

119. 0 

1. 05 

1.79 



He 

137. 0 

1.20 

1.83 

5 

2.64 

He 

20.0 

0. 17 

.78 



He 

41.6 

.30 

1.23 



He 

61.6 

.40 

1.36 



He 

77.8 

.45 

1.44 


i Average Mach number over center half inch of stream. 
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not show this general trend. In these nozzles the pressure 
first decreased, then increased, as the probe was moved from 
the nozzle wall to the center of the stream. Figures 1 1 and 
12 show this effect. In general, the static -pressure variation 
across the stream was not large. The average change over 
the center l / 2 inch of the stream amounted to 2 percent, 
while the maximum deviation was only 5 percent. 

The impact-pressure surveys of all the nozzles had one 
common characteristic; that is, the impact pressure increased 
from a minimum value at the edge of the nozzle wall to a 
maximum value near the center of the stream. The area of 
the stream which had a constant, or nearly constant, impact 
pressure varied somewhat from nozzle to nozzle, and also 
to some extent with pressure level. The average variation 
of impact pressure over the center ){ inch of the stream for 
all nozzles tested was 7 percent with a maximum variation 



Figure 11. — ' Variation of impact aud static pressure with nozzle radius for uozzlc number 1 

usiug nitrogeD gas. 


in some cases of 14 percent. Figures 11 and 12 show the 
variation of impact pressure across the stream for one 
nozzle using two gases, nitrogen and helium. 

The Mach number of a gas stream can be computed from 
measured values of the stagnation and static pressure of the 
stream, provided that the flow is isentropic between the 
points of pressure measurement. However, if the static 
and impact pressures are measured at the same point in the 
stream, no assumption regarding the nature of the flow need 
be made in order to compute the local stream Mach number. 
For an inviscid flow, the Mach number can be calculated 
from Rayleigh’s equation, the use of which involves the 
assumption of a normal shock wave ahead of the impact- 
pressure tube and isentropic compression of the gas in the 
subsonic flow behind the shock wave to the stagnation point 
on the impact-pressure tube. In a low-density flow, how- 
ever, viscous effects are present in the flow about a stagna- 
tion region, thus invalidating the assumption of isentropic 
compression behind the normal shock wave. Chambre 
(reference 14) computed the viscous correction to be applied 
to spherical-head impact tubes in a low-density flow, arid 
figure 13 shows the Mach number distribution of nozzle 1 
as computed by the Rayleigh equation (inviscid) and as 
computed by Chambre’s equation. The Mach numbers cal- 
culated from the stagnation and local static pressures were, 



Figure 12.— Variation of impact aud static pressure with nozzle radius for nozzle number 1 

using helium gas. 
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Figure 13 .— Comparison of Mach number distribution in nozzle number 1 as calculated 
with and without viscous correction, test gas, nitrogen. 


in all cases, higher than those computed using local values 
of impact and static pressure, thus indicating that the flow 
through the nozzle was nonisentropic due to viscous effects. 
Chainbre’s equation was used to compute Mach number 
from measured static and impact pressures for the tests 
described herein. 

The variation of impact and static pressure across the 
stream resulted in a stream in which the Mach number was 
not constant. Figures 13 and 14 show this effect and are 
generally representative of all nozzles tested. If only the 
center )z inch of the nozzles is considered, the average 
variation in Mach number amounted to 5 percent; however, 
at certain pressure levels in some of the nozzles the variation 
was as much as 14 percent. 

In order to determine if shock waves were present in the 
center % inch of the gas stream, the nitrogen-glow-discharge 
technique of flow visualization was used to examine the 
flow qualitatively. This technique was originally described 
in reference 15. The method used to generate the glow 
discharge was to pass the test gas through a screen which 
was connected to a pulsed direct-current power supply. 
The screen was located at the entrance plane to the nozzle 
and was pulsed to a negative peak potential of 1000 volts. 
The pulsation frequency was 1000 cycles per second. 

Figure 15 shows the gas stream from nozzle 1 at a static 
pressure of 114 microns. A 5° half- angle wedge was placed 
in the stream to generate a relatively weak shock wave 
which could be used as a yardstick to evaluate the strength 
of the nozzle shocks, if present. No shock waves were 



Figure 14.— Comparison of Mach number distribution in nozzle number 1 as calculated with 
and without viscous correction, test gas, helium. 

observed in the gas stream at this pressure level if the wedge 
was removed. As the static-pressure level was lowered, the 
shock wave from the wedge became too indistinct to be 
photographed. Observation of the nitrogen glow did not 
indicate the presence of strong shock waves in the nozzles. 

DRAG AND HEAT TRANSFER TESTS 

The test procedure used in the drag and heat-transfer 
experiments was to calibrate the drag balance with the 
tunnel shut down and then to set tunnel conditions corre- 
sponding to those obtained during tunnel-calibration tests by 
matching the stagnation and test chamber static pressures. 
As extremely close control over tunnel conditions was diffi- 
cult to attain, the pressure level during certain test runs 



Figure 15. — Shock wave from 5° half-angle wedge at pressure of 114 microns as seen with 
nitrogen glow discharge. 
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varied approximately 5 percent from those obtained during 
the tunnel-calibration tests. Examination of the results of 
the tunnel-calibration tests showed that Mach number varia- 
tion is small for small changes in the pressure level; therefore, 
the Mach number obtained during nozzle-calibration tests 
was applied to the heat-transfer and drag tests without 
correction. The Mach number used in plotting the data 
(which was converted numerically to molecular-speed ratio 
by means of the relation s = M ^yj2 ) was the average value 
obtained from a graphical integration of the Mach number 
distribution across the center l % 2 inch of the stream. 

I After stable operating conditions had been attained, the 
1 drag force, model temperatures, and other pertinent data 
I were recorded. Table II lists the recorded data. It will 
be noted that, in the column listing the recorded drag force, 
f two values occur on some of the runs. The different values 
j were due to separate drag measurements during the course 
of one run or resulted from small shifts in the calibration of 
the drag balance as determined before and after each run. 
In the figures which are later presented, showing drag coeffi- 
cients, the points correspond to those calculated using the 
I average of the two values listed in table II. 

To determine if the model surface conditions were changing 
with time, the model drag and temperature were measured 
at K-hour intervals over a period of 6 hours for one particular 
test condition. Figures 16 and 17 show the results of these 
tests and it can be seen that the drag coefficient and tempera- 
ture ratio remained essentially constant, thus indicating that 
the accommodation coefficient and tunnel conditions were 
constant over this period of time. 

The reproducibility of the data was checked by repeating 
the above-described run on separate da}^s. The reproduc- 
ibility was found to be within the accuracy of the drag and 
temperature measurements described above. 


RESULTS AND DISCUSSION 


A DISCUSSION OF THE ASSUMPTIONS ANI) RESULTS OF THE THEORETICAL 

ANALYSIS 

An examination of the heat-transfer equations (14) and 
(15) for a monatomic and a diatomic gas, respectively, reveals 
several interesting features. If an adiabatic case is consid- 
ered, that is, if both internal-energy input to the cylinder 
and radiant-energy exchange are considered absent, it can 
be seen that the ratio of cylinder temperature to free-stream 
static temperature is a unique function of the molecular- 



Figure Hi.— Variation of measured drag coefficient with time. 
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Figure 17.— Variation of the ratio of model to free-stream temperature, T e /T with time. 


speed ratio. The values of the temperature ratio computed 
from equations (14) and (15) are shown in figure 18. For 
comparison, curves of the temperature ratio for continuum 
flow are also shown. These latter curves were computed for 
the case where the total stream energy is recovered — in 
other words, the temperature recovery factor was considered 
to be unity. Several remarkable conclusions can be drawn 
from an examination of this figure. It can be seen that an 
insulated cylinder in a free-molecule-flow field will attain a 
higher temperature at a given value of molecular-speed ratio 
(or Mach number) than will be obtained in a continuum 
flow field. Since the continuum curves also represent the 
temperature ratio of a gas that has been accelerated adia- 
batieally from rest to any value of molecular speed ratio, it 
may be inferred that an insulated cylinder located in a low- 
density wind tunnel under free-molecule-flow conditions will 
attain a temperature which exceeds the total temperature of 
the gas stream. This result, of course, is in direct contrast 
to the corresponding phenomenon which occurs under con- 



Figure 18.— Variation of the ratio of cylinder temperature to free-stream statietemperature, 
with molecular speed ratio for an adiabatic cylinder in free-molecule flow. 
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tinuum-flow conditions where an insulated body can have, 
at most, a temperature equal to the stream stagnation temper- 
ature and normally does not even attain this temperature due 
to outward heat flow in the boundary layer. The anomaly 
can be explained, however, by a consideration, in the case 
of free-molecule flow, of the magnitudes of the incident and 
re-emitted molecular energy. The incident molecular energy 
is computed using the total velocity resulting from combina- 
tion of the stream mass velocity and the random thermal 
velocities. When the total velocity term is squared there 
results a term, which is twice the scalar product of the vector 
mass velocity and the vector thermal velocity. This scalar 
product term effectively increases the apparent internal 
energy of the gas from the continuum value of 3/2 kT per 
molecule to a value varying from 2 kT to 5/2 kT depending 
upon the speed and orientation of the body. The apparent 
internal energy becomes 5/2 kT for high speeds and body 
angles of attack greater than zero. This is just equal to the 
internal plus potential energy per molecule of a cube of 
continuum gas. Therefore, the energy incident on the body 
in free-molecule flow becomes equal to that of a continuum 
for large values of the molecular speed ratio for surfaces 
inclined at angles of attack greater than zero. The molec- 
ular energy which is re-emitted from the surface is assumed 
to be equal to the energy of a stream issuing effusively from 
a Maxwellian gas, in equilibrium at some yet unspecified 
temperature, into a perfect vacuum. For the case of an 
insulated body, the temperature of this gas is that of the 
body. The energy of the re-emitted stream calculated in 
this manner is equal to 2 kT w per molecule. The correspond- 
ing energy for a continuum gas at the same temperature is 
5/2 kT w per molecule if both thermal and potential energies 
are considered. For a given re-emitted stream temperature, 
a smaller amount of energy per molecule is transported from 
a body for the case of free-molecule flow than is transported 
in the case of continuum flow. It is clear that, if the same 
total amount of energy, namely the incident energy, is to be 
removed in each case, the effusive stream temperature and, 
hence, the body temperature for an insulated body must be 
higher for free-molecule flow than for continuum flow. 

It is seen from figure 18, that, theoretically, a higher cylin- 
der temperature is predicted for the case of a monatomic gas 
than for a diatomic gas. The explanation lies in the addi- 
tional energy which the diatomic gas is able to remove from 
the cylinder by virtue of the internal energy component due 
to molecular rotation. Of course, the impinging diatomic 
gas transmits more energy to the cylinder due to the rota- 
tional energy component; however, the rotational component 
is a smaller fraction of the total impinging energy than of 
the total re-emitted energy. 

An examination of the equations for the drag coefficient of 
an insulated cylinder (equations (17) and (18)) shows that 
the total drag coefficient, C D =C Di +C Dr , is a function only 
of the molecular-speed ratio. A plot of the total drag co- 
efficient as a function of the molecular-speed ratio for an 
insulated cylinder is shown in figure 19. It is notable that 
the drag coefficient approaches infinity as the molecular- 
speed ratio approaches zero. Hence, if the expression for 
total drag is written in the conventional drag-coefficient 



form, it can be shown from equations (17) and (18) that C D 
becomes inversely proportional to the speed ratio as $ 
approaches zero. It is also notable that the drag coefficient 
approaches a constant value as the speed ratio becomes large. 
This means, simply, that the effect of the random thermal 
motions becomes negligible as the speed of the cylinder 
relative to the gas becomes high. 

HEAT-TRANSFER TEST RESULTS 

The salient jpoint of the analysis concerning heat transfer 
to a cylinder in a free-molecule field was the prediction that 
an insulated model would attain a temperature in excess of 
the tunnel stagnation temperature. This result was observed 
in all runs. The maximum rise of the model center-point 
temperature above the tunnel stagnation temperature was 
65° F, using nitrogen as the test gas, and 147° F, using 
helium as the test gas. 

The results of the heat-transfer tests are shown in figures 
20 and 21 in which the ratio of the cylinder-center tempera- 
ture to the free-stream static temperature is plotted as a 
function of the molecular-speed ratio for the nitrogen and 
helium runs, respectively. The stream-static temperature 
was computed from the adiabatic relationship, using the 
measured value of Mach number (or molecular-speed ratio) 
corresponding to the nozzle pressure ratio as determined 
from the nozzle-calibration tests. The temperature-ratio 
curve for continuum flow assuming a temperature- 
recovery factor of unity is shown for comparison. It will 
be noted that in all cases the temperature ratio exceeded the 
maximum which could be attained under continuum-flow 
conditions. It can be seen that, in general, the test points 
indicate that the cylinder did not attain the full temperature 
rise that was predicted by the free-molecule analysis. This 
result may be ascribed to heat losses occurring from radiation 
and from conduction along the wire to the end supports. 


A COMPARISON OF THEORY AND EXPERIMENT FOR HIGH-SPEED FREE-MOLECULE FLOW 
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The fact that a large temperature gradient existed along the 
length of the wire due to end losses is clearly shown by the 
cylinder-end temperatures which are listed, together with the 
cylinder center-point temperature, in table II. 

DRAG TEST RESULTS 

The measured drag force on the cylinder was reduced to 
conventional dimensionless drag-coefficient form by means 
of the equation 



of drag coefficient. One of the basic assumptions of the 
drag analysis is that the drag coefficients are independent 
of Reynolds number or Knudsen number. In contrast, the 
continuum- and slip-flow regimes are characterized by the 
dependence of drag coefficients both upon Reynolds number 
and Mach number (cf., reference 16). Although it was not 


C D 


D 


D 


yLrpM 2 2 Lrps 2 


( 20 ) 


? 




The data points obtained from equation (20) are shown 
in figure 22. For comparison, the theoretical curves of drag 
coefficient as a function of molecular-speed ratio for an 
insulated cylinder, taken from figure 19, are shown as a 
solid line. It can be seen that agreement between theory 
and experiment is good considering the limitations of the 
theory and the difficulties attendant with the measurement 
of drag forces from 1 to 20 milligrams. The experimental 
values for the drag coefficient tend to fall below the theo- 
retical values at molecular-speed ratios above about 1.0 
and to be higher than the theoretical values at speed ratios 
below 1.0. The former characteristic may be due to a tend- 
ency toward specular reflection at the higher speed ratios 
while the latter characteristic, at subsonic values of speed 
ratio, may have been caused by interference with the flow 
caused by the shields at the ends of the model. 

The major premise of the existence of free-molecule flow 
was borne out by several factors in addition to the good 
agreement between the calculated and experimental values 



Figure 22.— Variation of total drag coefficient with molecular speed ratio. 
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experimentally possible to vary the Knudsen number inde- 
pendently of the speed ratio for one test gas, the use of 
both nitrogen and helium gave overlapping values of speed 
ratio with an approximate 2:1 variation of Knudsen number. 
The variation of Knudsen number covered by the tests with 
molecular-speed ratio is shown in figure 23. No trend is 
noticeable in the drag coefficients obtained in each case; 
consequently, it can be concluded that free-molecule flow 
existed. 

A comparison between measured drag coefficient and the 
drag coefficient calculated by methods as given by Lamb 
(reference 17) illustrates the inapplicability of continuum 
theory in the free-molecule regime. A sample calculation 
made by Lamb’s continuum method for conditions of molec- 
ular-speed ratio equal to 0.5 and a Knudsen number of 185 
(based on cylinder radius) gave a drag coefficient of 390. 
The corresponding experimental value was 14.6. 

CONCLUSIONS 

The following conclusions may be drawn from a compari- 
son of the test results with free-molecule-flow theory shown 
in this report: 

1. The salient point of the heat-transfer analysis, which 
was the prediction that an insulated cylinder would attain 
a temperature higher than the stagnation temperature of 
the stream, was confirmed by the test results. 

2. Drag coefficients, calculated from free-molecule-flow 
theory, agreed well with the measured drag coefficients. 
As predicted theoretically, the measured drag coefficients 
were independent of the Reynolds and Knudsen numbers and 
depended primarily upon the molecular-speed ratio. 

Amcs Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 12 , 1950. 


APPENDIX A 

THE MOMENTUM AND ENERGY EXCHANGE PROCESSES ACCOMPANYING MOLECULAR EMISSION FROM A SURFACE 


For the case of diffuse scattering of a stream of molecules 
from a surface, all directional history of the impinging stream 
is erased upon impact with the surface. The impinging 
stream is assumed to be scattered in a random manner which 
will be described in detail later. The energy and momentum 
carried by the scattered stream depends upon the velocity 
of emission which can be related to the energy (or tempera- 
ture) level that the stream has. The effective temperature 
or energy of the scattered stream depends upon the efficiency 
of the energy-exchange process that occurs between the solid 
surface and the impinging stream. This efficiency can be 
described by introduction of the accommodation coefficient 
which is defined as 


a= 


Ei-E r 

E t — E w 


(Al) 


where the symbols have the following definitions: 

Ei rate of incident molecular energy 
E T rate of re-emitted molecular energy 
E w rate of re-emitted molecular energy that would be carried 
by the scattered stream if it were a stream issuing from 
a gas in equilibrium at the surface temperature, T w 
It is assumed that the emitted stream has a Maxwellian 
distribution of speed corresponding to a gas in equilibrium 
at an unspecified temperature, T t . This fictitious gas has, 
also, an unspecified number of molecules per unit volume, 
N r . Now, as shown in reference 8, page 62, the number of 
molecules scattered from unit surface area per second that 
have speeds in the range dv and (hat are moving in a direction 
that makes an angle lying in the range d9 with the normal 
to the surface is given by 


dn r =2irN r A T v z e~ 0r2r2 sin 9 cos 9dvd9 (A2) 

Then, the normal component of force imparted to the 
surface by this group of molecules is 

dG T ~2TNrA T mv i e~ Pr2r2 sin 9 cos 2 9dvd9dA (A3) 

and the translational energy carried by this group is 

d E r = TrN r A r 7 nv 5 e~ p r 2v2 sin 9 cos 9dvd9d A (A4) 

By integration of these expressions with respect to 9 and v, 
the following values are obtained for the total momentum 
and energy carried by the scattered stream from dA: 

d G r = 2 xN T A,mdA C v'e-^dv f "'^in 0 cos 2 ndQ=^ 0 dA 
Jo Jo 4/3/ 

(A5) 

dE r =TrN T A r mdA f v 5 e~^^ 2v2 dv f sin 9 cos 9d9= dA 
J o Jo 2 a / 7 r/3 r 3 

(A6) 


In order to obtain a value for N r in terms of known variables, 
the number of molecules scattered from the surface is equated 
to the number of molecules striking the surface. The num- 
ber of molecules per second scattered from unit area is 

n T =2TN r A r [ v*e~Pr 2 » 2 dv f sin 9 cos 9d9=— ^l - (A7) 

Jo Jo 2 V 7T/3 r 

and the number of molecules striking unit area per second is 
shown in reference 4 to be given by 


n= 


N x 

2^(3 


(A8) 


The value of N r is obtained by equating n r and n from equa- 
tions (A7) and (A8), thus, 


Nr=N X 


fir 

fi 


(A9) 


By substitution of this value of N r in the momentum and 
energy equations, the following is obtained: 


da ’~ N m. iA 

(AlO) 

d E, — A L mx dA 

2 V 7T fifir 2 

(All) 


An identical analysis yields a value for E w of 
d E w = Nn {* - A = 2nkT w dA 


(A 12) 


For the case of a diatomic gas where the expression for the 
re-emission energy per molecule contains an additional 
term having the value kT Wy which accounts for the rotational 
energy per molecule, the total re-emission energy has the 
value 3nkT w dA. 

Equation (AlO) describes the normal re-emission pressure 
on the front side of a unit inclined area. In order to obtain 
the net drag force (or pressure) in the direction of the mass 
velocity Z7, the pressure on the rear side of the inclined area 
must be taken into account. The computation of the pres- 
sure on the rear side of the area is made in a similar manner 
to the computation of the pressure on the front side of the 
area and leads to an expression 


dG/- 


Nm x ' 


4 fifir 


dA 


(A 13) 


Then, the net drag force due to re-emission pressure in the 
direction of Z7, the mass velocity, is 


dG-dG V 


Nm sin 6 

4 fifir ~ 


(x-x')dA 


(A 14) 
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Replacing x and x' by their values, equation (A14) becomes 

(A 15) 




2Mr 

This expression can be put into drag-coefficient form by use 
of the following substitutions: Nm = p , s/(3=U. Equation 
(A15) then becomes 


dG-dG r '=(^P) 


pU 2 \ V^r sin 2 6 


dA 


(A 16) 


where s r is the molecular-speed ratio of the re-emitted stream. 
Then the re-emission drag coefficient, based on the total 
area of one side of the plane area, is 




N 7 r sin 2 6 


(A 17) 


In order to obtain an expression for s r which contains 
measurable or calculable physical quantities it is necessary 
to write an equation expressing an energy balance on the 
body in question. For the case of an insulated cylinder, if 
radiant energy exchange is neglected, only molecular ener- 
gies are involved and the incident molecular energy can be 
equated to the re-emitted molecular energy. To write an 
energy balance equation for a transverse cylinder, it is 
assumed that the energy balance is maintained for the cylin- 
der as a whole and not for individual elements of the cylinder. 
This means, essentially, that the cylinder has a uniform 
temperature, circumferentially. The total incident molecu- 
lar energy on the cylinder is, for a monatomic gas, 

[IT" x(s2 +^ 0+ V(s 2 +iw] (A 18) 

The total molecular energy emitted from the cylinder surface 
is 


The values of the integrals are shown in appendix B to be 


9) 


j: 

j; 


■'/* z t +z, . 

xde= — s — ? +yirs 


0 

x/2 


'■ f de = 


Zj\ + Z 2 

= — o ■' 7rs 


x 

x/2 


/; 


X^dd= Z x -\-— Z 2 -\~x^Tr S 


xV'd0=Z,+| Z 2 — |vxs 


When Ei and E r are equated from equations (Al 8 ) and 
(A19) and the values of the integrals are substituted, the 
following result is obtained: 

i L[.>(Z, + Z ! )+ 4 -^±^]=^(z, + Z,)(A20) 


Therefore, 


Pr = P 


2 (Z, + Z 2 ) 


Z,(«* + 2 )+ Z 2 (s 2 +|) 


(A 21 ) 


A similar analysis for the case of a diatomic gas where the 
incident molecular energy is 

a -sM ,,,+ 1 x, ' #+<s ' + 1 >f ’ x ' ,l,+ 

/’x/2 fir/2 “1 

J x +dd+ j x’+'de J (A 22 ) 

and the re-emission energy is 
gives a value of of 


(A23) 


3(Zi + Z 2 ) 


"Y Z 1 ($ 2 +3) + Z 2 ^-s* 2 +^ 


(A24) 


Now, the total re-emission drag force acting on dA has been 
shown (cf., equation (A 16)) to be 

To obtain the total re-emission drag on a cylinder, equation 
(A16) must be integrated over the cylinder surface, thus, 

Or- G /={-— ) (2 rL) £ sill 2 ddd=^~-) (2rL) ^ 

(A25) 

Then, the re-emission drag coefficient based on the pro- 
jected cylinder area is 




*3/2 

4s, 


(A26) 


The final equation for the re-emission drag coefficient for an 
insulated cylinder in a monatomic gas stream is obtained by 
using the relation s,= [//3, and taking 13, from equations 
(A 21 ) and (A24). 


3/2 I Zi(s 2 + 2 )-\-Z 2 (s 2 +A 

— 2( z;w — (A27 > 

and, for an insulated cylinder in a diatomic gas stream, as 


C D — 


4 s 


'Zj(s 2 +3)-f-Z 2 ^ 


3(Zj +Z 2 ) 


(A28) 



APPENDIX B 

HEAT TRANSFER TO A TRANSVERSE CIRCULAR CYLINDER IN FREE-MOLECULE FLOW 


It has been shown (see Analysis) that an energy bal- 
ance on a cylinder located in a stream of monatomic gas 
flowing with a velocity U perpendicular to the cylinder axis 
is expressed by 

2 «rL[ ffn ( - d 2 -+#r)d<M- (~ + rkT) <**]- 

4 otrLk T c ( nd0+ ^^~n , dd^+2irrL[Q — ea(T c i —T, i )] = 0 

(Bl) 

Equation (Bl) can be simplified by dividing through by 
2 ttjL and replacing n and n f by their values from equations 
(2) and (4), thus, 


Nv tt 


4tt 3/2 


“(m[/ 2 -4A*r c )( f’ r \de+ \l‘\'dey 


akTNv , 

2ir 3/2 


( f;\*de + T,*)+Q= 0 (B2) 

The integrals to be evaluated are 

r*i 2 

j. 

C t /2 

X 

) o 

X '/2 


(B3) 

(B4) 

(B5) 


(B6) 


If the value of x in equation (B3) is replaced by its defi- 
nition (from list of symbols), equation (B3) becomes 

sin sin 8erj{s sin B)]dO (B7) 

The first term of this integral can be evaluated by setting 
s = 2x and cos 20= z 
dz 


Then, dO = 


2 sin 26 


and the limits become 


2 = 1 when 0=0 


z = — 1 when 0=— 

Zd 


Thus we can write 


6" 52sin2 ^/0=V —= dz 

Jo 2 .)_! ^/j _ 0 2 


(B8) 


Now, by definition (cf., equation (11), p. 46, reference 18), 
the modified Bessel function of order n can be written as 

; - w — vfrm’-Hi m (D’L <B9) 


Then 


and since 


i 0 (x)=- r ~ 1 — f —iZ 

vVr(i/2) J -i vi — 


r (|) = ’ v/7r 

the first term of integral (B3) becomes 


J e-» 2 / 2 / 


•G) 


The second term (cf., equation (B7)) is 


/- r x/2 

\tt s I sin 0d0=\irs 

Jo 


(BIO; 


(Bll) 


The third term of equation (B7) is 

- C T/2 

y ts I sin 6 erf ( s sin 0) dO 

Jo 

which can be integrated by parts. Let 
erf (s sin 0)=u 


and 

Then 

and 


sin 0d0=dv 


2s 

du = —= 6“* 2 sin 2 0 cos Odd 

yV 


v = —cos 0 

Then, the integral becomes 


r >/ 2 r t/2 

A 7 rs —cos 0 erf ( s sin 0) +2cs* 2 I cos 2 0e~ s2 s™ 2 °d0 

(B12) 

The bracketed term of equation (B12) is zero, leaving the 
second term yet to be evaluated. Make the substitutions 


cos 2 6= 


1 +cos u 


sin 2 0= J 


1 —cos u 


Then 


and, replacing limits 


s 2 =2x 




8=0 u=0 


1 9=~ 


U=TT 
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there results 

2s 2 f'^cos 2 9e- s2sin2e de= S ~ j^(l + cos 


Values of the dimensionless variables and Z 2 were com- 
puted from the tables of Bessel functions in reference 19. 
Then, equation (B21) is obtained by substitution of equations 
(B15) to (B20) into equation (B2). 


s 2 e x r* 

~~2~ Jo ' 


iU du+- 


*> p- x 


j: 


cos ue XC0SU du (Bid) 


Nv vi a 
4**'* 


(mU 2 -4kT e )(Z l + Z 2 )+ 


Integral (B13) becomes, when substituting 2 = cos u and 
changing limits, 


s 2 e z 


2 





dz 


ct k TNr 
2 id ' 2 


(2Z,+|z 2 )-e<r(7V-yV)+Q=0 


(B21) 


If mv m 2 /2=kT, kTN—p, and s=U/v m are substituted into 
equation (B21), the following equation results: 




(B14) 


which is obtained by integrating the second term of (B14) 
by parts and substituting equation (B9). Then, the final 
value for integral (B3) is 


r 


‘\de=le-*i*L 


°(0 


I — \ TTS 

+ + 


(BIS) 


Now, integrals (B4), (B5), and (Bb) contain only terms 
which can be evaluated by the methods shown above by 
which integral (B3) was evaluated. The results are 



x ' dO =^ e ~ s 2 l 2 Io ( j , )- e <2/2 [ /o (|) +/ ’(l)] 

(B16) 

| x^® =,re-s2/2 ^o ^J+2 V^ s + 

5 ,rs2r,2/2 [ 7 °(V) + 7i (^)] (B17) 

| xV , dfl=ire _,2/2 / 0 ^^—^y7rs+ 

^ S 2 e -^[/o(^)+/,@] (B18) 


I 11 order to eliminate an excessive amount of writing, new 
variables are defined as follows: 

Z^ve-WIo^) (Bl9) 

Z,= .«"e-*[/.(|')+ /,(£)] (B20) 


2 j'(Z,J-Z,)-[z,(«’+2)+Z,(s ! +|)]+ 

— [«<r(7V — T,*) — Q] = 0 (B22) 

J)Vm<x 

which is the final energy exchange equation for a monatomic 
gas. An identical treatment of equation (13) for a diatomic 
gas gives 


3 (Z, + Z 2 )-[z,(.y 2 +3)+Z 2 (W£)]+ 

— MTV- 77)- 0 = 0 (B23) 

pv m Oi 

The groups j (s) = Z x (.s 2 + 3) + Z 2 ^ s 2 + 4 Q, g (s)= 3 (Zj + Z 2 ), 

/°(s)=Z 1 (s 2 +2)+Z 2 and g°(s)=2 (Zj + Zo) are tab- 

ulated in the following table: 


s 

/(«) 

9 ( S ) 

f ° (*) 

9 ° (*) 

0.2 

9. 80271 

9. 612*34 

6. 59860 

6. 40823 

.4 

10. 94784 

10. 16410 

7. 55981 

6. 77607 

.6 

12. 89399 

11.04924 

9. 21092 

7.36616 

.8 

15.69746 

12. 22282 

11.62319 

8. 14855 

1.0 

19. 43609 

13. 63288 

14.89180 

9. 08859 

1.2 

24. 20781 

15. 22758 

19. 13197 

10. 15172 

1.4 

30. 12863 

16. 96068 

24. 47507 

11.30712 

1.6 

37. 32923 

18. 79408 

31.06454 

12.52939 

1.8 

45. 95230 

20. 69872 

39. 05273 

13.79915 

2.0 

56. 14932 

22. 65358 

48. 59807 

15. 10239 

2.2 

68. 07771 

24. 64392 

59. 86307 

16.42928 

2.4 

81.89947 

26. 65990 

73.01284 

17. 77327 

2.6 

97. 77902 

28. 69470 

88.21412 

19.12980 

2.8 

115.88287 

30. 74382 

105. 63493 

20. 49588 

3.0 

136. 37821 

32. 80398 

125. 44355 

21.86932 

3.2 

159. 43361 

34. 87300 

147. 80928 

23. 24867 

3.4 

185.21769 

36. 94914 

172. 90131 

24. 63276 

3.6 

213. >9927 

39. 03112 

200. 88890 

26. 02075 

3.8 

245. 64801 

41. 11800 

231.94201 

27.41200 

4.0 

280. 63340 

43. 20904 

266. 23039 

28. 80603 

4.2 

319. 02474 

45. 30354 

303. 92365 

30. 20236 

4.4 

360. 99204 

47. 401 C 8 

345. 19168 

31.60072 

4.6 

406. 70528 

49. 50126 

390. 20486 

33. 00084 

4.8 

456. 33364 

51.60366 

439. 13242 

34. 40244 

5.0 

510. 04756 

53. 70808 

492. 14487 

35. 80539 

5. 2 

568.01730 

55.81428 

549.41254 

37. 20952 

5. 4 

630.41182 

57. 92200 

611. 10449 

38. 61467 

5.6 

697. 40229 

60.03114 

677.39191 

40. 02076 
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APPENDIX C 

THE DRAG DUE TO IMPINGING MOMENTUM ON A TRANSVERSE CYLINDER IN FREE-MOLECULE FLOW 

The differential momentum imparted, by the impact of impinging molecules, to the front and rear sides of a body 
located in a free-molecule flow field is (cf., references 1, 2, and 3) 

dGi = ^ | j* J (lzC x +lyCy+l z c z )c x e-^ Cx ~ Ux ) 2+ ( c *'~ u ^ 2 +( c *~~ u *W dc x dc v dc 2 — 

( ( ( + l v c v + l z c z )c x e ^ Cx ~ u *y+^v- u v) 2 +(c*~ Uz) 2 ] dc x dc y dc x \dA 

In equation (Cl), l x , l v , / 2 are the cosines of the angle between 
the coordinate axes and the direction of the momentum 
force. Now let 6 be the angle of attack of the element, that 
is, the angle between the mass velocity U and the plane of 
the plate. Then the drag force is in the direction of U and 
the following relations can be derived: 


(Cl) 


U X =U sin 6 


U v = — U cos 6 


u z = 0 


Z x =sin 6 
ly= — cos 0 

l 2 = 0 


(C2) 


Then, integration of equation (Cl) and substitution of the 
relations of equation (C2) yield 

in _P#T 2 «~* 2sin2 « . 

dGt 2 Ly /2 s + 

2 sin 6 erf(s sin 6) (* + 2 ^)J^ (C3) 

It is remarkable that a purely analytic solution for this type 
flow leads to a simple equation of the type of (C3) in which 
the bracketed term is a dimensionless function of molecular 
speed ratio and which may be interpreted as an ordinary drag 
coefficient. 

In order to obtain the drag force on a transverse cylinder 
due to impinging momentum, equation (C3), which was 
derived for an elemental plane area, must be integrated over 
the surface of the cylinder. To this end, dA is written as 


dA = rLdd 


(C4) 


Then, for the cylinder, the impinging drag force is 

„.._ 0 T (pV 2 \ r 2 Posin’* ( 

Gt 2rL \ 2 )J 0 L *r 1/2 S + 

2 sin d erj{s sin 6) ^ 1 + 2^)1 ^ 


(C5) 


This integral can be evaluated by the methods shown in 
appendix B. The final value for the drag force due to 
impinging molecules is 


(pU^ 


bO 

V 

' * r°v 


<«> 

Thus, there results from equation (C6) an impinging drag 
coefficient based on the projected area of the cylinder, which 
has the value 


C k 


(C7) 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal 

X 

X 

Rolling 

L 

Y >Z 

Roll 


U 

V 

Lateral 

Y 

Y 

Pitching.. . 

M 

z — >x 

Pitch 

r 

Q 

v 

r 

a 

Normal 

z 

z 

Y awing 

N 

X >Y 

Yaw 

\b 

w 

r 








Y 



Absolute coefficients of moment 

C, ~q bS 

(rolling) 


r _M_ 

Cm ~qeS 

(pitching) 


o-JL 

n ~qbS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


D 

V 

P/D 

V' 

V. 

T 

Q 


4. PROPELLER SYMBOLS 


Diameter 
Geometric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient C T 


Torque, absolute coefficient C Q = — 


pn 2 D* 

_Q 


c, 

v 

n 

3> 


eYl 

1 Pn 2 


Power, absolute coefficient @p— 

Speed-power coefficient = ^ ^ 
Efficiency 

Revolutions per second, rps 
Effective helix angle=tan _1 ( 


pn 2 !? 

5. NUMERICAL RELATIONS 


F 

\2irrnJ 


\ 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower =0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m= 5,280 ft 
1 m=3.2S08 ft 
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